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ABSTRACT 

We present initial results of a 750 ks Chandra observation of the remnant of Kepler's supernova of AD 1604. 
The strength and prominence of iron emission, together with the absence of O-rich ejecta, demonstrate that 
Kepler resulted from a thermonuclear supernova, even though evidence for circumstellar interaction is also 
strong. We have analyzed spectra of over 100 small regions, and find that they fall into three classes. (1) The 
vast majority show Fe L emission between 0.7 and 1 keV and Si and S Ka emission; we associate these with 
shocked ejecta. A few of these are found at or beyond the mean blast wave radius. (2) A very few regions 
show solar O/Fe abundance rations; these we associate with shocked circumstellar medium (CSM). Otherwise 
O is scarce. (3) A few regions are dominated by continuum, probably synchrotron radiation. Finally, we find 
no central point source, with a limit ^100 times fainter than the central object in Cas A. The evidence that the 
blast wave is interacting with CSM may indicate a la explosion in a more massive progenitor. 

Subject headings: ISM: individual (G4.5+6.8 (SN 1604)) — supernova remnants — X-rays: ISM — super- 
novae: general 



1. INTRODUCTION 



tral analysis was performed with XSPEC v. 12 ( Arnaud 



Kepler's supernova of 1604 (Kep ler|1606| l was the most re- 
cent well-observed Galactic supernova. The supernova type 
has been controversial for decades; the Galactic latitude (b = 
6?8) implies a distance above the Galactic plane of 470 pc for 
a distance of 4 kpc (San krit et al.| |2005), arguing for a Type 
la (Population II) progenitor, but optical observations (Blair 
et al. 1991) show radiative shocks in dense knots, sugges- 
tive of circumstellar material (CSM) shed by the progenitor. 
Bandiera (1987) proposed that the progenitor was a massive 
runaway star with a strong wind; Kepler's bright N rim sug- 
gests a bow shock directed away from the Galactic plane. 

X-ray observations of young supernova remnants can give 
powerful constraints on progenitor models and explosion 
mechanisms; whether Kepler had a Type la or core-collapse 
origin, the detailed distribution of ejected material and the 
composition of CSM could be used to test SN models. We 
undertook a long observation of Kepler with the Chandra X- 
ray Observatory to study the composition and distribution of 
ejecta and CSM, to search for any compact object, and to 
search for convincing evidence of the SN type. 

2. OBSERVATIONS AND RESULTS 

We observed Kepler's SNR for 741 ks with the Chandra 
X-ray Observatory ACIS-S CCD camera as a Large Project 
in six installments between April and July 2006. In all six, 
Kepler was positioned (slightly differently) on the S3 back- 
illuminated chip. Data were processed using CIAO v3.4 and 
calibrated using CALDB v3.1.0. A large background region 
to the N of the remnant (covering most of the remaining area 
on the S3 chip) was used for almost all spectra except when 
local backgrounds were called for, as described below. Spec- 
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1996). 



We used the nonequilibrium-ionization (NEI) v2.0 thermal 
models, based on the APEC/APED spectral codes ( |Smith et] 
|al.|2001[) and augmen ted by addition of inner-shell processes 
( [Badenes et al.|2006| >. 

We obtained about 3 x 10 7 source counts, with fewer than 
3% background (though some of those may be dust-scattered 
source photons). S/N permitting, spectra were extracted be- 
tween 0.3 and 8 keV. Images were smoothe d using platel ets 
with a default smoothing parameter 7=1/2 ( Willett 2007). 

Figure [1] shows the full observation, with red for energies 
between 0.3 and 0.72 keV (including primarily O), green for 
0.72 - 1.7 keV (primarily Fe L-shell transitions, also includ- 
ing Ne and Mg Ka complexes), and blue for 1.7-8 keV 
(with a strong contribution from Si Ka). Immediately evident 
is the absence of dominantly soft-spectrum regions, and the 
hardness of emission at the periphery. A few regions do show 
reddish colors; representative spectra are discussed below. 

No central point source is obvious. The bright central knot 
is resolved, and several fainter central knots, though perhaps 
unresolved, also have spectra of shocked gas. A faint source 
could be concealed under bright extended emission. Our con- 
servative estimate assumes a high total background of 860 cts 
per ACIS S3 pixel, appropriate for the bright central band of 
emission. A 5a detection of a point source within the 80% en- 
circled power radius of 0."685 (6 px area) would require 400 
cts, or 6 x 10 -4 cts s -1 . The compact central object in Cas A 
( |Pavlov et al. ||2000| , if placed within Kepler and using the 
absorbing column density determined below, would produce 
0.18 cts s . We conclude that we would detect such a source 
even if its luminosity were ^100 times fainter. 

For spectral analysis, we fixed the absorbing column den- 
sity at Nh = 5.2 x 10 21 cm" 2 , a value we found by fitting line- 
less regions in the SE with an absorbed synchrotron model 
srcut using Wilms ~et al.|(|20 00) abundances. We used the 
three-color image of Fig. [1] to aid in identifying small, rela- 
tively homogeneous regions of contrasting spectral character. 
We identified and examined about 100 regions. The vast ma- 
jority show prominent lines of Fe L-shell transitions, He-like 





FIG. 1.— Merged image between 0.3 and 8 keV. Red: 0.3 - 0.72 keV; 
green, 0.72 — 1.7 keV; bl ue, 1.7 — 8 ke V. All three images were smoothed 
using platelet smoothing (Willett 2007 1. Image size is 4'.7 X 3 '9. Regions 
from which spectra of Figs.[2]and[3]were extracted are shown. 



Si, S, Ar, and Ca, and Fe Ka, and a handful had contrasting 
properties. We found that the spectra can be separated into 
three main classes. Over 90% of the regions showed Fe- and 
Si-rich spectra; we identify these with SN ejecta. A few re- 
gions showed spectra with O and Mg emission consistent with 
solar- or near-solar abundances; we identify these with CSM. 
Finally, a few regions showed continuum-dominated spectra. 

Examples of ejecta spectra (regions 1 - 3 on Fig. [T} are 
shown in Fig.[2]and the top spectrum (Fe L region, 6) in Fig. [3] 
In particular, Fig.[2]shows spectra of three outer knots, where 
we might expect to find the lightest nucleosynthetic products 
at the largest distances, specifically the O that should indicate 
a CC explosion. All of them, however, show strong Si, S, 
and Fe L emission; prominent Ar, Ca, and Fe Ka lines are 
also present in the bright N knot 1 . The western and southern 
knots 2 and 3 have comparable Si and S emission, but they 
differ in their Fe L-shell emission. 

Knot 3 is wel l modeled by a heavy-element plane shock 
( |Hendrick|2063| i with velocity v, of 1010 km s" 1 ([889, 1 120] 
km s -1 90% confidence interval), shock ionization age r = 
2.69(2.32,3.18) x 10 10 cm" 3 s, Si/Fe and S/Fe ratios (by num- 
ber) of 2.64 (2.35,2.92) and 2.24 (1.59,2.87), respectively. 
The knot is clearly composed of ejecta. A trace amount of 
Mg (Mg/Fe = 0.052 [0.003,0,10]) has been included, because 
it formally improves the fit, but its presence is much in doubt 
in view of deficiencies in Fe L atomic data. A synchrotron 
component (constrained by the absence of radio flux from the 
knot [DeLaney et al. 2002]) improves the fit quality. We as- 
sumed no collisionless heating at the shock; if such heating 
were present, v s could be substantially less. Si and S are the 
primary constituents of the southern knot, but a significant 
(1/4 by mass) amount of Fe is present. 

In the western knot 2, the shock speed is lower (v s = 
500[360,830] km s" 1 ) and the ionization age is longer 
(t = 4.8[2.7, 14] x 10 10 cm " 3 s) but most significant are 
much higher Si/Fe and S/Fe ratios, 11.6 (8.9,14.5) and 14.7 
(10.2,21.1), respectively. There is noticeably less Fe (6% by 
mass) in this knot; different Fe contents of knots 2 and 3 par- 
tially account for their different Fe L-shell emission strengths. 
Finally, in knot 1 the relative strengths of Si, S, and Fe L- 
shell emission are similar, so this is also a Si- and S-rich knot 



FIG. 2. — Spectra of three outer knots shown on Fig.[T] Top (red): Bright 
Si- and S-rich knot 1 in the N. Fe L shell, Fe Ka, Ar, and Ca emission are 
also strong. Middle (blue) and bottom (green): Knots in the S (3) and W (2), 
with plane shock model fits shown. 



containing Fe. (Because a simple plane shock model failed 
to provide a satisfactory fit at high energies, we have not at- 
tempted a detailed abundance analysis for knot 1). 

The striking lack of obvious O features in the vast major- 
ity of the regions we examined led us to search specifically 
for soft-spectrum regions that might be expected to have a 
different composition. We found a few (all showing a pink 
color in Fig. [TJ with discernible O features; two are shown in 
Fig. [3] They constitute our second class of spectral region. 
The central knot has optical and IR counterparts (Blair et al. 
199 1 1 |2007| l, implying that it is composed of ambient CSM. 
We obtained the spectrum shown in Fig. [3] after subtracting 
a local background from the bright E-W band nearby. In the 
plane-shock model fit shown in Fig. [3] we allowed the abun- 
dances of Ne, Mg, Si, and Fe to vary (with S and Ni tied to Si 
and Fe, respectively). The poor spectral resolution of CCD's 
and the uncertain calibration at low energies, combined with 
deficiencies in Fe L atomic data, do not allow us to deter- 
mine reliable absolute abundances with respect to H and He, 
or even a relative N abundance, but supersolar abundances for 
all elements except N are ruled out. We assumed near so- 
lar composition by setting the O abundance to its solar value, 
and determined relative abundances with respect to O. We as- 
sumed a supersolar (by a factor of 3; Blair et al. 2007 ) N abun- 



dance. We obtained a temperature kT = 0.64 keV (0.59,0.70), 
r = 2.7(1.8,3.6) x 10 11 cm" 3 s, and abundances with respect 
to solar of Ne = 0.64 (0.41, 0.91), Mg = 0.83 (0.61, 1.1), Si 
(tied to S) = 1.4(1.0, 1.9), and Fe (tied to Ni) = 1.1 (0.94, 
1.4). When N was allowed to vary, the best-fit N abundance 
is 3.2 (1.9, 4.6), where errors do not include systematic un- 
certainties particularly important at low photon energies. The 
model underestimates the continuum near the Si and S Ka 
lines, most likely because of the presence of multi temperature 
plasma in the central knot, so abundances of these elements 
may be somewhat overestimated. The central knot spectrum 
is consistent with the same N-enhanced near-solar abundance 
composition found from optical spectra of the CSM in Kepler. 
Fig.[3]also shows another soft-spectrum knot, 5; its spectrum 
is consistent with that of the central knot. 

Finally, a few regions, primarily in the SE, showed feature- 
less or almost featureless spectra (one example, region 4, is 
shown in Fig. [3]>. These are well-described by models of syn- 
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FIG. 3. — Spectra of four knots shown on Fig.[T] Top (blue): Fe L region 
6. No O or Mg features are evident. Second from top (green): Soft-spectrum 
knot 5. An inflection at 0.6 keV is O Lya, and a feature at Mg Ka is labeled. 
Third from top at 0.8 keV (red): bright central knot. Note O and Mg features. 
Bottom at 0.8 keV (magenta): synchrotron-dominated region 4 at rim in SE. 
Only a faint Si feature is apparent. 



chrotron emission extending from the radio (Reynolds et al. 
|2007| l; they form our third class of region, and will be de- 
scribed more fully in a future publication. 

3. KEPLER'S PROGENITOR TYPE 

Spectral regions 1, 2, 3, and 6 in our first, Si- and Fe-rich 
class are almost certainly ejecta, of a composition completely 
inconsistent with a CC event, but expected for a thermonu- 
clear explosion. In fact, they resemble Si- and Fe-rich com- 
pact ejecta knots in the v icinity of the eastern blast wave in 
Tycho ( Decourchelle et aL]|2001 



2. Scarcity of O. By contrast, there is little evidence for O 
in Kepler's SNR. It is seen in only a few regions with spectra 
quite different from the ejecta knots in having relatively weak 
Fe L emission. Most of these regions are coincident with op- 
tical knots which have been identified with CSM of roughly 
solar abundances except for elevated N. The X-ray knots do 
not require enhanced N abundance, but are consistent with the 
optical abundances. In particular, O/Fe is roughly solar. In SN 
ejecta of Type la, O/Fe is expected to be only 0.3-0.7 by num- 
ber, whereas in CC explosions, it is about 70 (e.g., Iwamoto et 
al. 1999). None of the spectra we examined had the high O/Fe 
ratios required for CC ejecta. (Several poorly-understo od lu- 
minous Type Ib/c SNe listed by Richardson et al. ( 2006 ) could 
be an exception, because their O/Fe ratios are much less than 
typical for CC explosions. Without detailed explosion models 
for these SNe confirming their CC origin, we cannot exclude 
such possible extreme CC events on the basis of the O/Fe ra- 
tio alone, but they are ruled out by arguments below.) We 
conclude that regions with detectable O are CSM. 

3. Central source. No central point source is seen, with a 
limit at least 100 times fainter than the source in Cas A. 

We find the X-ray evidence for Kepler's origin in a ther- 
monuclear event to be compelling. However, additional evi- 
dence can be found from optical observations. 

4. Balmer-dominated shocks. Blai r et aL] ( |1991| l 
discovered substantial emission from nonradiative, Balmer- 
dominated shocks. Such shocks demand the presence 
of partially neutral upstream material, ruling out most 
CC scenarios. A hot (i.e., compact) progenitor of a 
stripped CC Type Ib/c SN would have a strong ioniz- 
ing flux before the explosion, while a cooler star would 
be larger at the time of the explosion, resulting in a 
large ionizing flux at shock breakout. The presence of 
neutral H sets an upper limit to the progenitor radius of 
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the blast wave with such an advanced nuclear burning compo 
sition have not been found so far in any CC SNRs, even Cas 
A; except for the narrow jet, O is instead the most common el- 
ement in the outermost ejecta (Fes en et al.|2006| l. Only about 
5% of our regions fell into the second class of roughly solar 
O/Fe and perhaps elevated N, as expected for CSM. We can 
state with confidence that nowhere in Kepler are there spectral 
regions consistent with the expected CC O/Fe ratio of 70 (e.g., 
|Iwamoto et al.|1999") ; it is very difficult to locate O spectral 
features at all. Where strong Fe and Si indicate ejecta, O is ab- 
sent. However, the dense CSM shown in optical observations 
and now identified in a few regions in X-rays is not expected 
for a Type la progenitor. In view of the unusual nature of our 
result, we collect the evidence on Kepler's progenitor type. 

1. Integrated X-ray spectrum and ejecta. The integrated 
spectrum (which resembles region 6, Fig. [3] top), as well as 
most of our spectral regions, shows very strong Fe L, no ap- 
parent O or Mg, strong Si and often strong Fe K emission. 
These elements are all identified with ejecta. (O and Mg lines 
can be seen in the RGS spectrum; Ballet 2004.) Even knots 
1, 2, and 3 at the farthest periphery show evidence of strong 
nuclear processing. The spectra resemble the integrated spec- 
trum of the la remnant of Tycho's SNR (Badenes et al. 2006), 
but are totally unlike that of the young O-rich CC remnant, 
1E0 102-70.4, which is dominated by O and Ne emission lines 
(Sasaki et al. 2006). Previous work with ASCA and XMM- 
Newton (Kinugasa & Tsunemi 1999, Cassam-Chenai et al. 
2004), derived substantial masses of shocked Fe, 0.1-0.3 M Q 
( Kinugasa & Tsunemi|1999 1. 



( Chevalierj2 005), where 13.6a eV (a > 1) is the mean energy 



of an ionizing photon, and Mqsm and M e j are the sh ocked 
CSM and ejecta masses. Since M CSM is only ~ 1M Q (Blair 
et al.|2007 1, explosions of red supergiants with R* ~ 6OO/?0 



are ruled out, including not only Type IIP and Type IIL SNe 
with full and partially-stripped H envelopes, respectively, but 
also Type lib SNe like 1993 J with only a residual (O.2M ) H 
envelope left on its He core (Woo sley et al.|19 94j). While an 
exotic stripped CC event, resulting from a star with a residual 
hydrogen envelope of low (<C 0.2M Q ) mass but large radius, 
might avoid ionizing the CSM before or during the explosion, 
such events become even less probable than a Type la with 
CSM interaction. An additional problem with a stripped 
CC progenitor is a large expected mechanical luminosity 
of its stellar wind, since radio observations of circumstellar 
interactions in Type Ib/c SNe are consistent with WR winds 
(e.g., |Chevalier|2007] l. A powerful WR-like wind would have 
accelerated and swept away any dense CSM to distances of 
many pc from the SN, in obvious c onflict wit h observations. 

5. Light curve. Stephenson & Green (2002) summarize the 
historical data. Kepler's observed peak magnitude of m ~ -3 
gives My = — 18.8 (using the observed Eg-v = 0.9 to obtain 
Ay = 2.8), brighter than most CC events but typical for SNe la 
(e.g., |Hamuy et al. |199 6). The light curve, though resembling 
a normal la, could also describe a SN IIL (but see above). SN 
1987A-like events (fairly compact blue stars with massive H 
envelopes) produce light curves that are far too broad, and are 
subluminous as well. 
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4. DISCUSSION AND CONCLUSIONS 

We believe that X-ray and optical evidence points com- 
pellingly to the conclusion that Kepler resulted from a ther- 
monuclear event. However, our deep Chandra observations 
also confirm that Kepler's blast wave is encountering mate- 
rial with at least solar metallicity, coincident with optically 
emitting N-enhanced material - most likely circumstellar ma- 
terial lost by a fairly massive star, either the progenitor of 
the white dwarf or the companion, or both. Interaction with 
CSM is unexpected for Type la events. Some Type la SNe 
have shown signs of CS M interaction in the form of narrow 
H lines, like SN 2002 ic ( |Hamuy et al."p003] > and SN 2005gj 
( |Aldering et al.|2006] l, but the interpretation of these objects 
remains controversial (Benetti et al. 2006; but see Prieto et 
al. 2007). Patat et al. (2007) give spectroscopic evidence that 
SN 2006X is a normal Type la SN with CSM interaction. 
Evidence for CSM interaction in SNRs appears to be rela- 
tively more widespread, and in some cases can be associated 
with Type la events. Prominent Fe emission in DEM L238 
and DEM L249, two LMC SNRs, indicates a Type la origin. 
However, the brightness and location of that Fe, at small radii 
with a large ionization age, cannot be explained with stan- 
dard Type la explosions into a uniform ISM (Borkowski et al 
2006); interaction with a nonuniform CSM seems to be re- 
quired. Relatively bright central emission dominated by Fe is 
common in LMC SNRs ( |Nishiuchi et al. |2001 van der Hey- 
den et al.~y2 004), suggesting that this class of SNRs is not a 
negligible population. N103B fLewis et al.|2 003) also shows 
strong Fe and Si emission and evidence for CSM interaction. 
Kepler's SNR could just be the nearest example of a class of 
thermonuclear supernovae interacting with dense CSM from 
the progenitor system. 

The required CSM might indicate an origin in a younger, 
more massive population. Three recent studies of extra- 
galactic SNe la (jMannucci et al.|[2006l |Sullivan et aLl|2006] 
Aubourg et al. 2007) 1 propose the presence of a "prompt" pop- 
ulation (delay <5x 10 s yr) in addition to an "old" popula- 
tion (delay ~ Gyr) in later-type galaxies. If stars with initial 
masses less than about 8 M Q produce white dwarfs, the min- 
imum latency for a thermonuclear supernova would be only 
about 50 Myr. A substantial number of remnants of Type la 



supernovae should come from the "prompt" population. 

However, Kepler's origin remains problematic. The ele- 
vated CSM abundances rule out a very old populatio n, but the 
systemic velocity (> 200 km s'^ jBlair et al.|l99l"| seems to 
indicate either a halo population, inconsistent with the abun- 
dances, or a runaway, probably inconsistent with a binary 
system. (The Galactic bulge population would not extend to 
Kepler's Galactocentric distance of 4.5 kpc.) The total CSM 
mass of > 1M Q ( Blair et aL]|2007[ ) requires too large an ini- 
tial mass for an extreme halo population. Could Kepler have 
resulted from the thermonuclear explosion of a single star? 
Current conventional wisdom requires that all AGB stars eject 
most of their envelopes to leave a sub-Chandrasekhar white 
dwarf; if that ejection were occasionally impeded, a single- 
star Type la supernova might result (e.g., Tout 2005). If Ke- 
pler's progenitor were a single star liberated from the Galactic 
plane by the core-collapse explosion of its companion about 
3 Myr ago, its maximum lifetime would only be of order 50 
Myr, also arguing for a more massive progenitor. 

Our conclusion that Kepler resulted from a Type la event 
with CSM interaction has important implications. Models of 
SN la progenitor systems make widely varying predictions 
for the nature of any CSM (e.g., |Badenes et al.|200"7] >; these 
predictions may now be testable. Similarly, SN la explosion 
models can now be tested in detail against the abundances 
and spatial distribution of elements in Kepler, and can be con- 
trasted with Tycho's SNR, apparently the result of a conven- 
tional SN la encountering only low-density, more or less uni- 
form ISM (Bad enes et al.|2 006). The existence of a more mas- 
sive subclass of SNe la with properties that may be systemati- 
cally different has obvious important implications for cosmo- 
logical use of SNe la, as it predicts a redshift-dependent mix 
of the two types. Since the use of SNe la for cosmology re- 
quires correction of peak magnitudes from light-curve shapes 
with a relation calibrated locally, the application of that same 
relation to more distant SNe la may raise concerns. 
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